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Abstract 

To explore the kinetic adsorption under continuous and nonequilibrium states, an integration of 
continuous measurement and adsorption platform kinetics method was proposed, which was 
initially called the ICM-AP kinetics method, and a corresponding kinetic adsorption experimental 
method was developed. Adsorption experiments of europium (Eu) on Ca-bentonite, Na-bentonite, 
and the D231 cation exchange resin were performed using the ICM-AP kinetics method and 
continuous measurements. Because the kinetic experimental results observed in this study were 
different from those of traditional batch adsorption data, pseudo-first-order or pseudo-second- 
order kinetic models were unsuitable for fitting the experimental data. Hence, a liquid membrane 
diffusion (LMD) model was developed based on the assumption of simultaneous 
adsorption/desorption to discuss the mechanism of kinetic adsorption. The kinetic adsorption 
mechanism was also studied by using XPS. The results indicated that the proposed adsorption 
model can fit the experimental data more suitably, and the adsorption/desorption behaviors of Eu 
on bentonite and the D231 resin were simultaneously observed, suggesting that the adsorption 
kinetics of Eu(III[) was mainly dominated by hydrated Eu (III) ions on the liquid membrane. 
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1. Introduction 


Nuclear energy is receiving increasing attention from many countries as a clean and efficient source of 
energy. France, as a former nuclear power country, decided to restart nuclear power projects on a large scale and 
continue to maintain the status of a nuclear power country [1]. By 2021, China's nuclear power generation will 
account for 5.02% of its electricity generation for the first time, and this proportion will continue to steadily 
increase [2]. However, considering the large-scale use of nuclear energy, radioactive and highly toxic wastes have 
been released into the environment, causing an increased risk of damage to human health and ecosystems. 
Therefore, studies regarding the adsorption of radionuclides in ecological environments, including the 
development of models and experimental investigations, are critical for predicting the behavior of radionuclides 
released into the biological environment and for the safety assessment of nuclear waste disposal [3]. According to 


the Status and Trends in Spent Fuel and Radioactive Waste Management published by the IAEA [4], low-level 
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radioactive waste accounts for approximately 95% of the total radioactive waste, whereas high-level radioactive 
waste accounts for less than 1%. Meanwhile, reprocessing has been widely considered as an effective treatment for 


low-level radioactive waste, for which the separation of radionuclides from radioactive waste is critical [5-7]. 


Traditional adsorption is considered as an effective and widely adopted method for removing radioactive 
pollutants from wastewater [8-10]. Generally, research regarding the adsorption behavior of radionuclides 
primarily focuses on the following three main aspects: nuclide species, adsorbents, and thermodynamic/kinetics 
models characterizing the adsorption mechanism. The study of nuclide species, uranium, thorium, transuranic 
elements (Np, Am, Cm, and Pu), fission products (Cs, Cr, and Pm), and lanthanides (Eu, Sc, and Ce) has gained 
notable attention [11,12,13-16]. Regarding the adsorbents, several previous studies have focused on rock minerals 
(granite, attapulgite) [17,18], clay minerals (bentonite, and red or yellow soil) [19,20], metal-organic frameworks 
(MOFs) [21], and organic polymeric materials [22,23]; meanwhile, the extraction of uranium from seawater has 
received greater attention [24]. Investigations regarding the adsorption of Eu on clay minerals have practical 
reference values; for instance, owing to its strong adsorption capacity and weak water conductivity, bentonite has 
been considered a major buffer and backfilling material for the geological disposal of high-level radioactive waste 
[25]. Numerous thermodynamic and kinetic models have been developed to characterize and predict adsorption 
behaviors. For instance, a non-linear two-site Langmuir model was used to describe the kinetic adsorption. 
Previous studies [26,27] employed a two-site Langmuir model to investigate the adsorption mechanism of Cs and 
Se on mudrock. The non-linear heterogeneity-based isotherms have also been used to characterize the anaerobic 


and aerobic adsorption of Se and Cs on various clay rocks. 


In this study, among the various cations, europium (Eu) was selected to perform the relevant adsorption tests 
owing to its stable species, which was mainly the Eu** cation with a solution of pH<9. Varying the pH values 
would affect the concentration of H*, which would highly influence the adsorption capacity of Am(III) on Na- 
bentonite] [28]. Theoretically, the desorption fluctuation curves are more likely to be observed with a relatively 
higher temperature and initial concentration of Am(III). Hence, a series of tests for different pH values, 


temperatures, and initial concentrations were performed to investigate the fluctuation mechanism in this study. 


An adsorption model is typically constructed based on the assumption that adsorption reaches a state of 
equilibrium. In this case, the adsorption isotherm should be the same as the desorption isotherm because the theory 
of thermodynamic equilibrium assumes complete reversibility of the chemical reactions, thus a unique (C, Q) pair 
[25]. However, adsorption mechanisms are driven by various controlled reactions or physical phenomena that have 
the largest variability in the reaction times, ranging from a few seconds to several years [29]. The largest 
variability in the reaction times mainly depends on the activation energy E (kJ/mol), which is the difference in the 
activation energies of the adsorption E, and deadsorption E2. However, for kinetic adsorption experiments, 
previous studies [30,31] have focused only on adsorption or desorption kinetics. The time scale of the kinetic 
reaction in previous studies was relatively long, and the measurement methods (ICP-OES or ICP-MS) for the 
intermittent measurements of the samples had limitations. In this case, the kinetic curves always exhibited a single 
adsorption or desorption state. In fact, the adsorption and desorption kinetics should not be considered 
independently because they occur, simultaneously, where one of the two dominates the entire process. Therefore, 
the adsorption and desorption processes should not be treated separately. This study proposes a new kinetic 
experimental method corresponding to a continuous measurement method to satisfy specific experimental 
conditions. The integration of a continuous measurement and adsorption kinetics platform was designed. An 
exponential kinetics model of liquid membrane diffusion was developed based on these assumptions to simulate 
the kinetic experimental data obtained in this study and to calculate the activation energy E according to the 


Arrhenius law. It was assumed that the kinetic adsorption process can be divided into two steps as follows: 1) The 


hydrated Eu(II) ions in the solution diffuse to the surface of the liquid membrane of the adsorbent particles, while 
desorption also occurs in this area. 2) The adsorption of the hydrated Eu(II]) ions further occurs in the solid 
internal hole [20]. 


2. Theory 
2.1 A new liquid membrane diffusion model 


A new exponential kinetics model of liquid membrane diffusion was proposed based on the assumption of 
adsorption and desorption simultaneously occurring to describe the entire nonequilibrium adsorption kinetics of Eu 
on Ca-bentonite, Na-bentonite, and the D231 cation exchange resin. In this model, the adsorption of Eu on the 
adsorbent particles mainly occurred in the outer layer of the liquid membrane and internal holes. The entire 
adsorption and desorption process was divided into two steps. First, hydrated Eu(II) ions mainly diffused into the 
liquid membrane located at the outer layer of the bentonite or resin particles; the hydrated Eu(III) ions in the liquid 
membrane would then be released into the solution. Second, the hydrated Eu(II]) ions attached to the liquid 
membrane would further diffuse into the solid internal hole; simultaneously, the hydrated Eu(II) ions would then 
release from the solid internal hole into the liquid membrane (physisorption mainly controlled the entire adsorption 
process ). A schematic of this mechanism is shown in Figure 1. In the first step, the reaction equation for the solute 


diffusion into the liquid membrane can be expressed as follows: 


N+L>NL (1) 


Adsorption 


Liquid membrane TN +L=>NL 


N+S>NS 
Adsorption 
Step 2 ® 


a T 


© © Desorption @ 
aaka No NSSN+S 


Desorption 


Solid internal hole 


Figure 1: Schematic diagram of the kinetics of the adsorption mechanism. The complete kinetics adsorption 
process was presumably divided into the following two steps: 1) The nuclides diffuse into the liquid membrane, 
which can be expressed by N + L > NL; the nuclides also diffuse from the liquid membrane to the solution, 
which can be expressed by NL > N + L. 2) The adsorption of nuclides adsorbed in the liquid membrane on the 
solid internal hole will continue to simultaneously occur with desorption; these equations are expressed by NL + 


S> NS and NS=> N +S, respectively. 


In the aforementioned equations, N indicates the nuclide of the solute, L is the liquid membrane, and NL is the 
compound of nuclides that diffuses into the liquid membrane. On the other hand, the desorption process can be 
observed throughout the entire adsorption process; however, the desorption fluctuation curves observed were 
random, irregular, and hardly quantitative. Hence, only the empirical exponential rate of the adsorption reaction is 


given by Eq. (2): 


dC;/Co 
dt 


= ke” (2) 

where C, is the concentration of Eu in the liquid membrane (mg/L), Co is the initial concentration of the 
solution (mg/L), and t is the reaction time (s). kų and yı (1/T) are the coefficient and rate constant of the 
exponent equation, respectively. In the second step, the adsorbate (Eu) located at the liquid membrane diffuses into 
the solid internal hole of the bentonite and resin (the adsorption kinetic reactions are driven by the chemical 
surface complexation and caution exchange, respectively), for which the reaction equation can be expressed as 


follows: 
NL+S=WNS (3) 


Here, S is the solid internal hole and NS is the compound of nuclides adsorbed onto the solid internal hole. The 


reaction rate is given by Eq. (4): 


dC2/Co 


2? = Kae V2! (4) 


where Cy, is the concentration of Eu in the solid internal hole (mg/L); C1 and C3 are the initial concentrations 
of Eu in the solution (Co). Compared to the previous step, the nonequilibrium adsorption kinetics were relatively 
weak. Considering the entire adsorption kinetics process, both steps occurred simultaneously; however, the weight 
factors of the two steps in the entire process were different. If the experimental conditions were changed, the 


weight factors of the two steps would also change. Hence, the total reaction rate can be expressed by Eq. (5): 


dC; /Co 4 dCz/Co 


di dt =a XxX k,e %'+(1 — a) x k,e Y2" (5) 


Eq. (5) can be reshaped to Eq. (6) as follows: 


= d (%42) = [ax k;e™™1t+(1 — a) x k,e™”2t]dt (6) 


(0) 


The indefinite integration of the total reaction rate in Eq. (6) is given by Eq. (7): 
=> C,/Co=f[a x k,e7¥1'+(1 — a) x kpe ¥2"]dt (7) 


The indefinite result of the integral can be expressed by Eq. (8): 


= CyCo=ax(—) e-V1t4(1 — a) x> e7Y2t (8) 


npr À k k , 
The parameter a indicates the weight factor of the first step; T. and a. are constants, which are regarded as 
~¥1 ~¥2 


B and C. Reshaping this equation can be expressed as follows: 


c = a x BeV1' + (1 — a) x Ce-¥2" (9) 
0 


where C+ is the concentration of Eu in the solution at various reaction times (mg/L) and Cg is the initial 
concentration of Eu (mg/L). However, numerous unknown parameters significantly affect the accuracy of the 
fitting results. According to the indefinite integral derivation process, parameters “B” and “C” in Equation (9) 
should be constants, both of which can be considered as “1”, with parameters “k1” and “k2” equal to rl and r2, 
Ka? « 


respectively. Because parameters r1”, and “r2” were sufficient to accurately describe the experimental data, 


parameters B and C are not shown in Tables (1), (2), and (3). 


Nevertheless, the random, irregular, and hardly quantitative desorption fluctuation curves observed for CyCo are 
shown in Figures 6, 7, and 8, respectively; the desorption rate equations for the two steps in this study could not be 
determined. The Arrhenius law provides the relationship between the kinetics constant k and activation energy E 
(kJ/mol) [31], which can be used to calculate the standard energy of the two steps proposed in the exponential 


kinetics model of liquid membrane diffusion. 


k= Derr (10) 
Here, D indicates a dimensionless collision probability factor; R and T indicate the molar gas constant (8.314 
J/(K/mol)) and thermodynamic temperature (K), respectively. The “k” parameter in Equation (10) is based on the 
reaction constant. Using Equation (10), we calculated the activation energy E for steps 1 and 2. The “k” parameter 
in Equation (10) was directly obtained with the value of parameter “y1”, which is shown in Equation (8) and is not 
mathematically deduced. The AF calculated according to Equation (10) was only a relative value and provided a 
reference for the two steps. The calculated parameter of D was considered for the value of the weight factor a, 
which can approximate the collision probability instead of parameter D. The kinetics model of the liquid 
membrane proposed in this study is a non-linear equation; its parameters were directly fitted by MATLAB 


(version 2019b) instead of linearization, which was then estimated using a non-linear approach. 


2.2 Kinetics adsorption conditions 


To meet the specific experimental requirements (relatively short kinetics adsorption time within a maximum of 60 
min) of the refit thermodynamic adsorption platform, a new kinetics adsorption method is proposed. The main 
difference of the newly proposed kinetics method in traditional batch tests is the simultaneous adsorption of Eu on 
bentonite and, more recently, the Eu concentration observed by ultraviolet spectrophotometry. The kinetics 
adsorption times of the batch tests (the kinetics adsorption of the nuclides on bentonite generally occurs within 7 


days) were also greater than that of the ICM-AP kinetics method. 


Table 1. Various experimental conditions 


Initial Eu The volume of Ca- 
Expt. nos. Curves pH Temperature l l 
Concentration bentonite used (30g/L) 
(a) 4 25° 1.5 mg/L 0.15 mL 
Exp. 1 (b) 4 25° 1.5 mg/L 0.30 mL 
(c) 4 25° 1.5 mg/L 0.45 mL 
(a) 4 25° 1.5 mg/L 0.30 mL 
Exp.2 
(b) 5 25° 1.5 mg/L 0.30 mL 


(c) 6 25 1.5 mg/L 0.30 mL 
(a) 4 29” 1.5 mg/L 0.30 mL 
Exp.3 (b) 4 35° 1.5 mg/L 0.30 mL 
(c) 4 45° 1.5 mg/L 0.30 mL 
(a) 4 25° 1.0 mg/L 0.30 mL 
Exp.4 (b) 4 25° 1.5 mg/L 0.30 mL 
(c) 4 25° 2.0 mg/L 0.30 mL 


Tables 2. Detailed experimental conditions for the adsorption of Eu on the D231 resin (Figure 8) 


Initial Eu 
Expt. nos. Curves pH Temperature The volume of D231 (30¢/L) 
Concentration 

(a) 4 25" 1.50 mg/L 5.00 mL 
Exp.1 (b) 5 25° 1.50 mg/L 5.00 mL 

(c) 6 25" 1.50 mg/L 5.00 mL 

(a) 4 30° 1.50 mg/L 5.00 mL 
Exp.2 

(b) 5 30° 1.50 mg/L 5.00 mL 

(a) 4 40° 1.50 mg/L 5.00 mL 
Exp.3 

(b) 5 40° 1.50 mg/L 5.00 mL 


2.3 Integration of continuous measurement and adsorption platform kinetics method 


Generally, each CyCo curve was obtained approximately 2-3 times because if only obtained once, the fluctuation 
observed in the curves may have been caused by a signal interference. For instance, if air was sucked into a rubber 
tube, the adsorbance value would suddenly increase, thus affecting the results. Many key experimental conditions, 
such as the rotating speed of the peristaltic pump and the pH value of the mixed solution, influence the value of the 
absorbance, which further causes an apparent error in the value of C;/Co. Therefore, after a series of tests, the most 
suitable rotating speeds of the magnetic rotor and inlet peristaltic pump were maintained at 1500 and 20 rpm, 
respectively; the pH of the mixed solution was measured to be 2.7. A schematic of the experimental setup is shown 


in Figure 2. 
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Figure 2. Schematic of the integration of continuous measurements and kinetic adsorption. Notes: the solution 


flow direction is from left to right. 


The kinetic adsorption devices used in this study differed from traditional equipment, as well as the batch or 
column methods. As shown in Figure 2, the adsorption reactor was filled with a standard europium solution at 
concentrations of 1, 1.5, and 2 mg/L, and material powders of Ca-bentonite, Na-bentonite, and the D231 cation 
exchange resin with a solid/liquid ratio of 30 g/L. Note, the different solid/liquid ratios of the adsorbent will 
influence the experimental results in the kinetic adsorption process, whereas the solid/liquid ratio in this special 
condition hardly changes because the particle size of the material powders is sufficiently small to cause the entire 
solution to present a homogeneous state, indicating a proportional reduction of the solution and clay particles. To 
meet the requirements of continuous measurements, traditional visible spectrophotometry was performed in this 
study. When traditional spectrophotometry is used to test the concentration of the samples, many specific 
conditions must be satisfied, including the optimal pH of the samples and adequate time to achieve color-rendering 
stability. To meet these requirements, two electromagnetic heating stirrers (Figure 2) were used to adjust the 
temperature of the sample solution and to shorten the color rendering time. An integral kinetics experimental 
process can be divided into two major steps or parts (the adsorption reaction and real-time online detection and 
recording). As shown in Figure 2, the mixed adsorption solution in beaker 1 contained either Ca-bentonite, Na- 
bentonite, or the D231 cation exchange resin powders and a standard europium solution. A magnetic rotor was 
used to enable the adsorption reaction to fully react in the beaker, and a filter membrane with a diameter of 0.45 
um was connected to the nozzle of the injection rubber pipe to separate the sorbent powders from the europium 
solution. The port of the outlet rubber conduit was connected to a three-way valve to mix the arsenazo III solution 
containing acetic acid-sodium acetate and the Eu sample solution. Subsequently, the mixed solution was allowed 
to flow into a heating and cooling beaker. Finally, the mixed solution flowed into the modified cuvette to detect 
and record the absorbance of the mixed solution displayed in the computer program, and then flowed into the 
beaker of the waste liquid. Many key experimental conditions, such as the rotating speed of the peristaltic pump 
and pH value of the mixed solution, influence the value of the absorbance, which further causes an apparent error 
in the value of C:/Co. Therefore, after a series of tests, the most suitable rotating speeds of the magnetic rotor and 
inlet peristaltic pump were maintained at 1500 and 20 rpm, respectively, and the pH of the mixed solution 
measured was 2.7. Appropriate rotating speeds ensured that the adsorption reaction of europium on the adsorbents 
and color rendering occurred adequately. A rotation speed of 20 rpm for the peristaltic pump ensured an adequate 


adsorption reaction time. The initial Eu concentration was adjusted to maintain the absorbance observed in the SP- 


721E spectrophotometer within the range of the linear working area. Generally, the concentration of europium 


should not exceed 2 mg/L; the pH of the arsenazo III reagent was maintained at 2.7. 


3 Materials and method 


3.1 Materials 


All the raw materials and reagents used in the experiments were analytical grade without further purification 
except for the raw material of the D231 resin, which had a particle size of approximately 0.45 um and was 
ground and passed through ~100 mesh sieves prior to use. The bentonite sample was obtained from Gaomiaozi and 
converted into Na-bentonite and Ca-bentonite according to a previous study [31]. The Eu solution used in this 
study was purchased from Guobiao (Beijing) Testing & Certification Co, Ltd. The acetic acid-sodium acetate pH 
buffer solution with a value of 2.7 that was controlled by sodium hydroxide was purchased from Xilong Scientific 
Co, Ltd. The arsenazo III reagent was purchased from Xilong Scientific Co, Ltd., which had a concentration of 1.0 
g/L. An ultraviolet spectrophotometer (SP-721E), peristaltic pump (L100-1S-2), pH meter (PHS-3C), and X-ray 


photoelectron spectroscope were used for testing. 


3.2 Preparation of the integration of continuous measurement and adsorption kinetics platform 


Traditional adsorption kinetics research methods, such as batch or column experiment techniques, have been 
frequently used to study various kinetics-related questions; however, these methods have certain inherent flaws, 
such that the adsorption kinetics experimental data are intermittently measured. Because the measurement and 
adsorption reactions were not operational and simultaneously occurred, many studies [20,32] have suggested that 
the fluctuation of the adsorption capacity versus the reaction time observed in the adsorption kinetics curves was 
caused by measurement or artificial errors. Moreover, owing to the lack of integration of the adsorption and 
measurement equipment as well as the corresponding test methods, the irregular fluctuations observed in the 
kinetics curves have not been thoroughly investigated. Traditional kinetic adsorption methods cannot explain the 
irregular fluctuations observed in the early stages of kinetic adsorption (within one hour). The majority of Eu was 
adsorbed on the bentonite within one hour, which is indicated by the frequent fluctuations in the kinetics curves. 
Therefore, to accurately understand the adsorption kinetics, the integration of continuous measurements and 
adsorption/desorption kinetics was proposed; the traditional SP-721E spectrophotometer (purchased from INESA 
Analytical Instrument Co., Ltd.) used in this study was retrofitted to achieve a continuous injection and real-time 
online display of absorbance on a computer. The original cuvette was pierced at the bottom; a diagram of the 


physical device is shown in Figure 3. 


Figure 3. Physical image of the refitting cuvette 


For real-time online detection and recording, a visualization Chinese window program was developed, as shown in 
Figure 4. Using this software, experimental data can be collected from the aforementioned instruments; that is, the 


curve of the adsorption capacity versus time (calculated in terms of absorbance). 


Figure 4. Schematic of the visualization Chinese window program. 


3.3 Kinetic adsorption tests and modeling 


The kinetic adsorption method used in this study differs from that used in traditional batch experiments. Three 
types of adsorbent materials were prepared, and the operating concentrations used for this adsorption test were 
designed using a standard solution of Eu at concentrations of 1, 1.5, and 2.0 mg/L. The experiments were 
conducted at varying temperatures (25°, 35°, and 45°), and pH values (approximately 3.0, 4.0, 5.0, and 6.0) of the 
adsorption solution to investigate the different conditional influences on the results of the adsorption kinetics. The 
experimental data obtained using this specific kinetics method were simulated using a liquid-membrane diffusion 
exponent kinetic model to characterize the adsorption kinetics. The trust region reflective algorithm and the 
Levenberg-Marquardt option, which can potentially make the liquid membrane diffusion model effective, are vital 
tools for the reliable analysis of the experimental data obtained by a new adsorption kinetics method. The residual 
root mean square error (RMSE) was used to determine the goodness of fit and is defined by the following 


equation: 


RMSE = J (Qi — qi)’, (11) 


where Qj; is the experimental data obtained from the new kinetics adsorption method used in this study, q; is the 
fitting value corresponding with Q;, and m and v indicate the number of data points in the new kinetics 


experiments and the number of parameters within the liquid membrane diffusion exponent model, respectively. 


4 Results and discussion 


4.1 XPS study 


In this study, the kinetics of the adsorption of Eu(III) under different adsorption times was studied using XPS. As 
shown in Figure 5, the main peak position of Eu 3ds/2 varies at 15 and 30 mins, indicating different binding 
energies of the adsorbed Eu(III) at variable adsorption times. Because the adsorption reaction of the liquid 


membrane of the Ca-bentonite particles dominated the initial adsorption process, the main peak of Eu 3d5⁄2 at an 
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adsorption time of 15 min may be attributed to the hydration of Eu(III) on the liquid membrane. As the adsorption 
time increased, the hydrated Eu(III) on the liquid membrane further diffused to the surface, and surface 
complexation occurred inside the hole sites, increasing the binding energy of Eu 3ds5/2. The main peak positions of 
the Eu 3ds/2 kinetic adsorption at 15 and 30 min were close and fixed at 1130.8 eV and 1132.1 eV, respectively. 
This indicates that the Eu(II) adsorption on Ca-bentonite gradually stabilized. The binding energy of the hydrated 


Euli) on the liquid membrane was the lowest, and its combination was relatively unstable. 
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Figure 5. (1) XPS spectra of Eu 3ds5/2 on the Ca-bentonite samples: (a) kinetic reaction time of 15 min and (b) 30 
min. The samples were prepared as follows: 20 mg/L of CEu(III), pH=5, and the weight of the samples was 3 


mg. 


4.2 Analysis of the C:/Co curves for the adsorption kinetics of Eu on Na, Ca-bentonite, and the D231 resin 
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Figure 6. (C:/Co) curves for the experimental data of the kinetics tests and fitting curves of the exponential kinetics 

model of liquid membrane diffusion for the adsorption of Eu on Ca-bentonite. Figure (a), (b), (c), and (d) represent the 

different volumes of Ca-bentonite used, pH, temperature and initial concentration of Eu, respectively (Notes: The 


detailed conditions of the experiments are listed in Table 1). 


Prior to an explanation of the results, note that the fluctuation in the CyCo curves above is not attributed to the 
signal noise. 1) If the fluctuation in the curves is considered signal noise, it should be observed throughout the 
entire CyCo curves and at the beginning. 2) As shown in Figure 6 (b), the CyCo curves for the adsorption of Eu(IID) 
on Ca-bentonite as a function of pH (pH=5, 6) were smooth and no fluctuations were observed, which is explained 
in the following states. A previous study [16] suggests a similar fluctuation in the kinetic curves of the adsorption 


of Am(III) on Na-bentonite as a function of the contact time. 


The adsorption experiments were performed as a function of the volume of Ca-bentonite used, pH, temperature, 
and initial concentration of Eu, and each kinetic test lasted approximately one hour. As shown in Figure 6 (a), as 
the volume of Ca-bentonite used increased, the adsorption rate gradually improved; meanwhile, smaller irregular 
fluctuations were observed in the backend of the Ci/Co curves. This phenomenon is attributed to the physical- 
desorption behavior of Eu(III) on the liquid membrane of the surface of Ca-bentonite. Note, the entropy change of 
the hydrated Eu(III) ions diffused to the surface of bentonite is positive (AS>0). The dehydration step requires high 
energy, which is provided by the reaction temperature (25°); thus, the reaction rate apparently improves as the 
temperature increases. This is considered as step “1” of the liquid membrane diffusion model, when physical 
desorption nearly occurs; this is the rate-controlling step that occurs at the beginning, which is indicated by the 
C/Co curves (within 15 min). Based on the assumptions of the liquid membrane diffusion model, the equation of 
step 1 (NL => N + L) that is shown in Figure 1 can be used to represent this physisorption kinetics behavior. The 
C,/Co curves as a function of the pH (4, 5, and 6) are shown in Figure 6 (b). The descending behavior of the Ci/Co 
curves with pH values of 5 and 6 were smooth, which indicates that the entire reaction process was dominated by 
adsorption whereas the physical-desorption of Eu(II) on the liquid membrane of the surface of Ca-bentonite was 
hardly occurring. There are two main interpretations for this. The concentration of the ionized H” decreased, which 
further weakened the competitive adsorption between Eu(III) and H* (the surface charge of bentonite is negative 
and the competitive adsorption will influence the process of Step 2 (N + S > NS) shown in Figure 1). However, 
the chemical reaction process of the surface complexation was dominant at these pH values, for which the 


equations can be expressed as follows: 
= SOH + Eu*+ == SOEu** + H* 


= SOH + Eu?* + H,0 = SOEuOH* + 2H* 


= SOH + Eu3* + 2H,0 == SOEu(OH)3 + 3H* 


The Eu adsorption process on bentonite does not change the physical structure of bentonite, rather only forms the 
complex compounds of SOEu**, SOEuOH*, and SOEu(OH)». As indicate in the chemical equations above, the 
equilibrium of the reaction moves forward. These two main factors ultimately led to a steady decline in the pH of 
the CyCo curves. When the kinetic adsorption tests were performed at different temperatures, as shown in Figure 6 
(c), the discrimination of these curves was subtle; a few irregular fluctuations in the backend of the curves were 
also observed. This phenomenon clearly indicates the physical-desorption of Eu(III) on the liquid membrane of the 
Ca-bentonite surface. As indicated by the Ci/Co curves in Figure 6. (d), (b), and (c), there were hardly any 
fluctuations, which is in agreement with the adsorption expectations of the physical-desorption of Eu(III) on the 


liquid membrane of the surface of Ca-bentonite owing to the excessive Eu(III) in the solution medium. 


Meanwhile, upon comparing the curves of the adsorption of Eu on Ca-bentonite, the CyCo curves for the 
adsorption of Eu on Na-bentonite were similar, as shown in Figures 7 (a), (b), (c), and (d), and only the adsorption 


rates were slightly different. 
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Figure 7. Experimental data curves (CyCo) of the kinetics tests and fitting curve of the exponential kinetics model 
of liquid membrane diffusion for the adsorption of Eu on Na-bentonite. Figure (a), (b), (c), and (d) present the 
different volumes of Ca-bentonite used, pH values, temperatures, and initial concentrations of Eu, respectively 


(Notes: the detailed conditions of the experiments are shown in Table 1.) 


The process of Eu adsorption on the D231 cation resin did not change the physical structure of the D231 resin; a 


cation exchange with the functional group of the D231 resin occurred. Irregular fluctuations at the backend of the 


C,/Co curves were also observed; however, the observed times were difficult to determine. These observations 


indicate that desorption occurs during the adsorption process, but do not confirm the precise time. However, as 


shown in Figure 8 (a), (b), and (c), the irregular fluctuations in the C:/Co curves for the adsorption of Eu on the 


D231 cation exchange resin are clearly observed, indicating that the desorption was relatively stronger than the 


adsorption on bentonite, and its adsorption capacity was lower than that of bentonite. The kinetic adsorption of Eu 


on D231 was completely driven by the chemical mechanism of the cation exchange with a relatively lower 


selectivity. Its adsorption combining energy was also relatively weak, which leads to a reversible reaction easily 


occurring. In contrast, the fluctuations observed in Figures 6 and 7 were relatively weak because bentonite has a 


higher adsorption capacity and selectivity for Eu and demonstrates a relatively weaker desorption phenomenon. 


(a) Exp.1 
1.O~ 
h, fa} pH=4 
Ei E =-=- (D) pH=5 
\ 7 ' (e) pH=é 
og LN hy —— The Fitting Curve 
Ma a E a (b) 
oe ii AF Err T ee 
LPs) ea te 
| mir (a) 
044 om I 
g2 = : z ; = ; = z : ; = ; 
[i] 400 Logn | 200 2000 7500 3000 
Time(s) 
(c) Exp.3 
1.0 
" (a) T=40°, pH=4 
ie -=+ {b} T=40°, pH=5 
—— The Fitting Curve 
of4 A 
cro 64 
šj r (a) 
Üd J n a Ricci 
=(b) 
g2 : ; s ; : ; s i = ; : ; s | 
o S00 1000 1300 2000 2500 3000 3500 
Time(s) 


(b) Exp.2 
| 0 | 
b, (a) T=30°, pH=4 
Dy --++(b) T=30°, pH=5 
| oan —— The Fitting Curve 
0.8 
0.64 = 2 (a) 
U | a i 
Fi 44 
(b) 
ies ri A La 
0.2 mente | 
0.0 . r n r : r : r : ' : r : | 
t s00 1000 1500 2000 2500 3000 3500 


Time(s) 


Figure 8. Curves of the experimental data (CyCo) of the kinetics tests and the fitting curve of the kinetics model 


of the liquid membrane diffusion exponent for the adsorption of Eu on the D231 cation exchange resin (Notes: 


The detailed conditions of the experiments are shown in Table 2). 


4.3 Analysis of the parameters of the exponential kinetics model of liquid membrane diffusion 


Table 3. Parameters of the exponential kinetics model of liquid membrane diffusion for the simulation of 


experimental data of the adsorption of Eu on Ca-bentonite 


a ¥1(1/T) 
The volume 0.15ml 0.4 6.89x10° 
used of Ca- 0.30ml 0.57 1.05x107 
bentonite 0.45ml 0.87 1.21x10? 


¥2(1/T) RMSE AE, (kJ/mol) AE, 
6.20x10°  1.45x10% 10.1 22,7 
1.26x10* 8.50x10° 9.89 20.2 
2.21x10* 6.68x10° 10.6 15.8 


a ¥i(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE> 
The pH of 
l 3 0.57 1.05x107 1.26x104 8.50x105 9.89 20.2 
mixed 
l 4 0.72 1.01x107 2.24x104 6.55x105 10.6 17.6 
adsorption 
l 5 0.79 4.66x10° 1.01x104 5.55x104 12.7 18.9 
solution 
6 0.71 7.56x103 7.60x104 2.09x104 11.3 14.7 
The temp of a ¥i(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
mixed 25° 0.57 1.05x107 1.26x104* 8.50x10° 9.89 20.2 
adsorption 35° 0.80 3.33x10° =1.23x10* 6.38107 13.6 18.3 
solution 45° 0.67 9.13x10° 2.31x10* 2.07x107% 10.7 18.0 
e a ¥i(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
e 
1.0ppm 0.57 1.05x107 1.26x104 8.50x105 9.89 20.2 
concentration 
E 1.5ppm 0.47 8.21x10° 5.56x10° 3.79x105 10.1 22.6 
of Eu 
2.0ppm 0.30 8.26x10° 2.49x10°  2.07x10°> 8.91 25.4 


Notes: Here, the parameters AF,and AE, (kJ/mol) indicate the activation energy of the adsorption of Eu on the 
liquid membrane and solid internal hole, respectively (details regarding the experimental conditions above are 


shown in Table 1). 


The RMSE values for the simulation of the exponential model of liquid membrane diffusion shown in Tables 3, 4, 
and 5 were considerably small, with values ranging between 1.45 10*—2.07x 10> (at least four or five orders of 
magnitude). These results indicate that the adsorption model proposed in this study is effective. The key 
parameters of the liquid membrane diffusion model and the calculated activation energies are listed in Tables 3, 4, 
and 5. As shown in the Table 3, as the volume of Ca-bentonite used increased, the parameter a representing a 
weight factor for Step 1 gradually increased from 0.40 to 0.87, and the values of the rate constant y4 were larger 
than yz (1.21107~6.89x 103>1.26x107~6.20x 10>). Moreover, the activation energy of Step 1 (AE,, the 
activation energy of the adsorption of Eu on the liquid membrane) was lower than AE, (the activation energy of 
the adsorption of Eu on the solid internal hole); a lower activation energy resulted in a faster and easier reaction 
rate. As the volume of Ca-bentonite used increased, the values of AE, hardly changed (9.89~10.6 kJ/mol)); 
however, the values of AE, gradually reduced (22.7—15.8 kJ/mol). These results can be attributed to the fact that 
the kinetic adsorption tests in this study were conducted in a nonequilibrium and unsaturated state. The total 
volume of the liquid membrane of Ca-bentonite gradually increased, which further led to more hydrated Eu(II) 
ions diffusing into the liquid membrane and, consequently an improvement in the value of the weight factor (a). 
Meanwhile, as the volume of Ca-bentonite increased, the surface complexation sites (= SOH) also rapidly 
increased, resulting in an easy performance of the complexation reaction. Therefore, the value of the activation 
energy AE, that was calculated according to the Arrhenius law shown in Table 3 gradually reduced (22.7—15.8 
kJ/mol). When the pH of the mixed adsorption solution progressively increased, the value of a increased and 
remained at approximately 0.7, indicating that the hydrated Eu(II) ions diffused into the liquid membrane (Step 1) 
and dominated the entire adsorption kinetics process. The rate constant y4 was larger than y2 (4.6610 
3~1.05107>1.01104~7.60 10%), and the activation energy AE, was larger than AE, (14.7~20.2>9.89~12.7 
kJ/mol). This also indicates that Step 1 occurred more easily and rapidly compared to Step 2. As the temperature 
increased, parameters a and AE, first increased and then decreased. Combined with the Ci/Co curves at a 
temperature of 35°, Step 1 for the hydrated Eu(IID) ion diffused into the liquid membrane (a=0.8) and dominated 
the entire kinetics process. When the temperature increased, the calculated AE, value gradually decreased from 
20.2 to 18 kJ/mol, and Step 2 easily occurred. The adsorption tests under different initial Eu concentrations were 


conducted; the trend for parameters a, AE,,and AE, were different from the other conditions. As the initial Eu 


concentration increased, the total concentration (Ci+Cz2) of adsorbed Eu on the liquid membrane and solid internal 
holes increased, whereas the volume of the liquid membrane did not change. However, the relative proportion 
concentration (C1) would decrease because Step 2 for the adsorption of Eu on the = SOH sites would further 
improve (the improvement of concentration (C2) was larger than C1), which caused the reaction rate to decrease 
(the gradual increase of AE (20.2~25.4 kJ/mol)). The key parameters for the adsorption of Eu on Na-bentonite 
are also listed in Table 4, which indicates that the adsorption kinetics process based on the assumptions of the 


proposed adsorption model theory is similar to the adsorption of Eu on Ca-bentonite. 


Table 4. Parameters of the exponential kinetics model of liquid membrane diffusion for the simulation of 


experimental data of the adsorption of Eu on Na-bentonite 


a ¥i(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
The volume 0.15ml 0.27 5.43x10° 7.14«10° = 2.80x«10> 9.71 22.9 
used of Na- 0.30ml 0.54 5.13x10° 1.77x104 5.16x«10° 11.0 20.0 
bentonite 0.45ml 0.77 4.19x10°  1.58xK107 8.77x105 12.0 17.6 
a yı(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AEF, 
The pH of 
3 0.54 5.13x10° 1.77x104 =2.80x«10> 11.0 20 
mixed 
4 0.44 7.31x10% 1.50107 3.49107 10.2 20.4 
adsorption 
5 0.62 7.92x10° 2.36x10* 4.19107 10.8 18.3 
solution 
6 0.63 6.16x10% 1.88x10* 1.46107 11.5 18.8 
The temp of a yı(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
mixed 25° 0.54 5.13x10° 1.77x104 2.80x10° 11.1 20.0 
adsorption 3a 0.62 7.18x10° 3.30x104 5.86x10°> 11.1 17.4 
solution 45° 0.63 5.68x10° 7.54x10° 1.30107 11.7 21.0 
a yı(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AEF, 
e 
l 1.0ppm 0.54 5.13x10° 1.77x104 2.80x10° 11.0 20.0 
concentration 
E 1.5ppm 0.45 9.58x10° 2.95x10° 2.68x105 9.55 24.4 
of Eu 
2.0ppm 0.30 1.03x10% 2.53x10° 3.02x10°> 8.33 25.4 


Notes: Here, the parameters AF,and AF, (kJ/mol) indicate the activation energy of the adsorption of Eu on the 
liquid membrane and solid internal hole, respectively (details regarding the experimental conditions above are 


shown in Table 1). 


The calculated reaction rate of Step 1 (y1) was always larger (4.19 10°~5.43 x 10°>1.58x 107~7.14 10-5) than 
Y2. This indicates that the front-end (rapidly decreasing) and rear-end (slowly decreasing) of the curves were 
owing to the adsorption of Eu on the liquid membrane and solid internal hole, respectively. The parameters for the 
adsorption of Eu onto the D231 cation-exchange resin are listed in Table 5. For the series of pH, the a parameters 
gradually increased (0.23—0.45) and the adsorption capacity of the D231 resin improved according to the Ci/Co 
curves shown in Figure 8. This is because the adsorption mechanism of Eu on the D231 resin was a cation 
exchange, which mainly depended on the pH of the solution. As the pH increases, the Eu(II) on the liquid 
membrane further supplies the demand of Step 2. The calculated reaction rates of Step 1 (y1) were always larger 
(5.13X103~7.92x 107>1.50x107~2.36x 107) than those of Step 2 (yz). The calculated weight factor a and the 
activation energy of the two steps were hardly influenced by the temperature. Comparing the different adsorption 
materials, bentonite and the D231 resin, the parameters of a for the bentonite were larger than those of the D231 


resin. Considering parameter of a for the D231 resin, the weight factor of a for the adsorption of Eu on 


bentonite was further influenced by the different experimental conditions. However, according to the calculated 
parameters of a shown in Table 5, the adsorption of Eu on the D231 resin was mainly controlled by Step 2 


(chemical behavior: cation exchange). 


Table 5. The parameters of the liquid membrane diffusion exponent kinetic model for the simulation of 


experimental data of the adsorption of Eu on D231 caution exchange resin. 


The pH of a y¥i(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
mixed 4 0.23 5.26x10°  1.57x10% 8.62 10> 9.36 21.1 
adsorption 5 0.30 2.03x107 7.29x10° 4.43x10° 12.4 22.7 
solution 6 0.45 2.36x10° 8.73x10°  5.63x10° 13.0 21.7 
The temp of a yı(1/T) ¥2(1/T) RMSE AE, (kJ/mol) AE, 
mixed 23” 0.23 5.26x10° 1.57x104 8.62x105 9.36 21.1 
adsorption 35° 0.29 2.00x10° 9.95x10°  5.08x10°> 12.3 22.0 
solution 45° 0.27 3.81x10° 1.95x104  4.74«10°> 10.5 20.4 


Notes: Here, the parameters AF,and AE, (kJ/mol) indicate the activation energy of the adsorption of Eu on the 
liquid membrane and solid internal hole, respectively (details regarding the experimental conditions above are 


shown in Table 2). 
5 Conclusions 


In this study, the desorption reaction was observed in the C:/Co curves when kinetic adsorption tests were 
conducted under different conditions. A new automated adsorption platform for the continuous measurement and 
its corresponding kinetic adsorption experimental method were proposed to conduct kinetic adsorption 
experiments of Eu on Ca-bentonite, Na-bentonite, and the D231 cation exchange resin. A new exponential kinetics 
model of liquid membrane diffusion was used to simulate the experimental data obtained to evaluate the 
characteristics of the adsorption kinetics of Eu. The model used for the adsorption kinetics in this study was 
consistent with the experimental data. According to the CyCo curves of Eu on Ca-bentonite, Na-bentonite, and the 
D231 resin, irregular fluctuations in the back-end of the curves were clearly observed, indicating the simultaneous 
occurrence of adsorption and physical desorption. The desorption of Eu on the D231 resin was more significant 
than that on the bentonite, which may be attributed to the lower selectivity of Eu on the D231 resin functional 
group. The values of a shown in Tables 3 and 4 were larger than those in Table 5, which indicates that the weight 
for Step 1 in the process of Eu adsorption on bentonite was stronger than that of D231; this implies that the surface 
charge of bentonite was negative, which improved the process of Eu adhering to the surface. Meanwhile, for the 
D231 caution resin, the surfaces of the resin particles were electrically neutral (the weight factor of step 1 was 
lower than that of bentonite), and the adsorption of Eu on the solid internal hole was carefully exchanged in an 
equal-charge manner. The calculated reaction rate of Step 1 (y1) was always larger than yz, which also reflects 
that the front-end (rapidly decreased) and rear-end (slowly decreased) of the curves were indicative of the 
adsorption of Eu on the liquid membrane and the solid internal hole, respectively; these results were in agreement 
with the CyCo curves and reflected the reaction activation energy AE. The XPS results demonstrated that the 
Eu(iII) adsorption on Ca-bentonite gradually stabilized. The binding energy of the hydrated Eu(III) on the liquid 
membrane was the lowest, and its combination was relatively unstable. A previous study [20] also indicated 


similar conclusions, which supports the results of this study. 


Based on the aforementioned analysis, the adsorption models proposed in this study can be integrated into 
chemical and physical adsorption models. The adsorption mechanism is shown in Figure 9, where the kinetic 
curves observed in this study were not smooth, but fluctuated, and curve (a) was integrated with curves (b) and (c). 
Curve (b) represents the chemical adsorption occurring on the surface and internal hole of the bentonite /D23 1 
resin, and curve (c) represents the physical adsorption occurring in the liquid membrane. Curve (b) represents the 
chemical adsorption types, indicating that the chemical bonds were stable and not immediately desorbed, whereas 
curve (c) represents the process by which the hydrated Eu ions diffuse into the liquid membrane, which is 
relatively unstable. Therefore, the amplitudes of the fluctuations in the Ci/Co curves observed in this study were 


controlled by the weights of the two steps. 


The whole adsorption curve observed were 
nat smooth but fluctuation and the (a) curve 
were Integrated with the (b), (c) curve. 


~ ‘The chemical adsorption occurred in the 


F surface and internal hole of materials. 
yp (C) 
Phe physical adsorption occurred in the 


liquid membrane of materials. 


Adsorption quantity 


Adsorption time 


Figure 9. Schematic of the adsorption types for the models proposed in this study. 


6 Acknowledgements 


This project was financially supported by the Natural Science Foundation of the Jiangxi Province, China 
(20202BABL203004), Opening Project of the State Key Laboratory of Nuclear Resources and Environment (East 
China University of Technology) (2022NRE23), and Opening Project of the Jiangxi Province Key Laboratory of 
Polymer Micro/Nano Manufacturing and Devices (PMND202101). 


7 Author contributions 

All authors contributed to the study conception and design. Material preparation, data collection and analysis were 
performed by Haoqi Yu, Tao Yu and Jianhua Ye. The first draft of the manuscript was written by Haoqi Yu and all 
authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. 


8 Reference 


[1] H. Wu, C.Y, Xu, X.H, Liu et al., Radioactive Waste Disposal-Key to Sustainable Development of Nuclear 
Energy. Nucl. S. 12, 155-159 (2013). doi:10.16432/j.cnki.1672-5360.2013.s1.028 


[10] 


[11] 


[12] 


J. Birkholzer, J. Houseworth, C.F. Tsang., Geologic Disposal of High-Level Radioactive Waste: Status, Key 
Issues, and Trends Annu. Rev. Environ. Resour. 37, 79—106 (2012). doi:https://org/10.1146/annurev- 
environ-090611-143314 

V.L. Loon, M. Glaus, C. Ferry et al., Studying radionuclides migration on different scales: the 
complementary roles of laboratory and in situ experiments. 446-483 (2012). 
doi:https://org/10.1533/9780857097194.2.446 

IAEA.Status and Trends in Spent Fuel and Radioactive Waste Management 
https://www.iaea.org/publications/1 1173/status-and-trends-in-spent-fuel-and-radioactive-waste-management 
M.F. Hamza, J.C. Roux, E. Guibal., Uranium and europium adsorption on amidoxime-functionalized 
magnetic chitosan microparticles. Chem. Eng. J. 344, 124-137 
(2018).doi.https://org/10.1016/⁄j.cej.2018.03.029 

S. Virtanen, S. Meriläinen, M. Eibl et al., Adsorption competition and kinetics of trivalent cations (Eu, Y and 
Cm) on corundum (a-A1203): a batch adsorption and TRLFS study. Appl. Geochem. 92, 71-81 
(2018).doi:https://org/10.1016/j.apgeochem.2018.02.011 

C.F. Wu, Y.W. Cai, L. Xu et al., Macroscopic and spectral exploration on the removal performance of 
pristine and phytic acid-decorated titanate nanotubes towards Eu(III). J. Environ. Chem. Eng. 6, 842—848 
(2018). doi:10.1016/j.molliq.2018.02.110 

R.S. Abolfazl, H.B. Ahmad, S.H. Hosseini et al., Kinetic, equilibrium and thermodynamic studies on 
adsorption of uranium and thorium from aqueous solutions by a selective impregnated resin containing 
carminic acid. J. Hazard. Mater. 286, 152-163 (2015). doi:https://org/10.1016/j.jhazmat.2014.12.047 

S.A.M. Idris, Adadsorption, kinetic and thermodynamic studies for manganese extraction from aqueous 
medium using mesoporous silica. J. Colloid. Interface. Sci. 440, 84—90 (2015). 
doi:https://org/10.1016/j.jcis.2014.10.022 

W.Q. Hao, B. Di, Q. Chen et al., Development of the Wade equation for the description of elution peak 
profile by using second-order rate equations for the adsorption kinetics. Chem. Eng. Res. Des. 96, 15-22 
(2015). doi:https://org/10.1016/}.cherd.2015.02.002 

Y.B. Cai, Y.T. Mi, H. Zhang, Kinetic modeling of antimony(III) oxidation and adsorption in soils. J. Hazard. 
Mater. 316, 102—109 (2016). doi:10.1016/).jhazmat.2016.05.027 

L. Xu, T. Zheng, S. Yang et al., “Uptake Mechanisms of Eu(III) on Hy droxyapatite; A Potential Permeable 
Reactive Barrier Backfill Material for Trapping Trivalent Minor Actinides”. Environ, Sci. Technol. 50, 
3852-3859 (2016). doi:10.102/acs.est.5b05932 

T. Yu, Y.Jun, Z.Y. Wang, Adsorption behaviour of Eu(III) on natural bamboo fibres: Effects of pH, humic 
acid, contact time, and temperature. Nucl. Sci. Tech. 31, 4 (2020). doi:10.1007/s41365-019-0710-3 

R.Q. Liu, X.P. Wang, S.Y. Ning et al., Adsorption behavior of *4!Am(IID) and Eu(IID) by silica/polymer- 
based isoHex-BTP adsorbent from nitric acid solution. Nucl. Sci. Tech. 26, S10301 (2015). 
doi:10.13538/j.1001-8042/nst.26.S10301 

F.L. Song, X.W. Yang, X.L. Li et al., The behaviour of cestum adsorption on zirconyl pyrophosphate. Nucl. 
Sci. Tech. 27, 60 (2016). doi:10.1007/s41365-016-0054-1 

B. Li, J.L. Liao, J.J. Wu et al., Removal of radioactive cesium from solution by zinc ferrocyanide. Nucl. Sci. 
Tech. 19, 88-92 (2008). doi:10.1016/S1001-8042(08)60029-9 

K. Fukushi, Y. Hasegawa, K. Maeda et al., Adsorption of Eu(III) on granite: EPMA, LA-ICP-MS, batch and 
modeling studies. Environ. Sci. Technol. 47, 12811-12818 (2013). doi1:org/10.1021/es402676n 

Z.S. Chen, J.T. He, L. Chen et al., Adsorption and deadsorption properties of Eu(IIJ) on attapulgite. J. 
Radioanal. Nucl. Chem. 307, 1093-1104 (2016). doi:10.1007/s10967-015-4252-9 


[19] 


[20] 


[21] 


[22] 


[23] 


[24] 


[25] 


[26] 


[27] 


[28] 


[29] 


[30] 


A. Chawla, M. Prasad, R. Goswami et al., Kinetic model for adsorption of divalent heavy metal ions on low 
cost minerals. Korean. J. Chem. Eng. 33, 649-656 (2016). doi:10.1007/s11814-015-0177-9 

T. Yu, Z.T. Xu, J.H. Ye, Adsorption kinetics of Eu(III) and Am(III) onto bentonite: analysis and application 
of the liquid membrane tidal diffusion model. J. Radioanal. Nucl. Chem. 319, 749-757 (2019). 
doi:https://org/10.1007/s10967-018-6386-z 

K.A, L. Chen, X. Dong et al., Efficient separation between trivalent americium and lanthanides enabled by a 
phenanthroline-based polymeric organic framework. Chin. Chem. Lett. 33, 3429-3434 (20272), 
doi:https://org/10.1016/j.cclet.2022.02.011 

A. Samokhvalov, Aluminum metal—organic frameworks for adsorption in solution: a review. Coord. Chem. 
Rev. 374, 236—253 (2018). doi.https://org/10.1016/j.ccr.2018.06.011 

M. Tuzen, T.A. Saleh, A. Sari et al., Interfacial polymerization of trimesoyl chloride with melamine and 
palygorskite for efficienturanium ions ultra-removal. Chem. Eng. Res. Des. 159, 353-361 (2020). 
doi:10.1016/j.cherd.2020.04.034 

J.H. Ye, T. Yu, Efcient and selective extraction of uranium from seawater based on a novel pulsed liquid 
chromatography radionuclide separation method. Nucl. Sci. Tech. 2, 13-27 (2023). 
doi:https://org/10.1007/s41365-023-01180-9 

D.G. Strawn, D.L. Sparks, Adsorption kinetics of trace elements in soils and soil materials. In: Selim, H.M., 
Iskandar,K.I. (Eds.), Fate and Transport of Heavy Metals in the Vadose Zone. Lewis Publishers, Boca Raton, 
FL, USA, pp.1—28. (1999) 

C.P. Lee, C.Y. Liu, M.C. Wu et al., Simulation of a 2-site Langmuir model for characterizing the adsorption 
capacity of Cs and Se in crushed mudrock under various ionic strength effects. J. Radioanal. Nucl. Chem. 296, 
1119-1125 (2013). do1:10.1007/s10967-012-2184-1 

H.Q. Yu, T. Yu, J.H. Ye et al., Simulation of non-linear heterogeneity-based isotherms for characterizing 
anaerobic and aerobic sorption of Se and Cs on various clay rocks. J. Radioanal. Nucl. Chem. 330, 1177-1189 
(2021). doi: https://org/10.1007/s10967-021-08012-3 

T. Yu, W.S. Wu, Q.H. Fan, Adsorption of Am(III) on Na-bentonite: Effect of pH, ionic strength, temperature 
and humic acid. Chinese. Chem. Lett. 23, 1189-1192. (2012). doi:10.1016/j.cclet.2012.07.011 

D.L. Sparks., Kinetics and mechanisms of soil chemical reactions. In: Sumner, M.E. (Ed.), Handbook of Soil 
Science.CRC Press, Boca-Raton, FL, USA, pp. B123—B167 (2000). 

E. Terte, G. Berge, E. Simoni, Europium retention onto clay minerals from 25 to 150 °C: experimental 
measurement, spectroscopic features and adsorption modelling, Geochim. Cosmochim. Acta. 70, 4563-4578 
(2006). doi:https://org/10.1016/j.gca.2006.06.1568 

M.B. McBride, (1994) Environmental Chemistry of Soils. Oxford University Press, New York, USA. (1994) 
T. Yu, W.S. Wu, Q.H Fan, Adsorption of Am(III) on Na-bentonite: Effect of pH, ionic strength, temperature 
and humic acid. Chinese. Chem. Lett. 23, 1189-1192 (2012). doi:10.1016/j.cclet.2012.07.011 


